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Magnetic domains at the surface of a ferroelectric monodomain BiFe03 single crystal have been 
imaged by hard X-ray magnetic scattering. Magnetic domains up to several hundred microns in size 
have been observed, corresponding to cycloidal modulations of the magnetization along the wave- 
vector fc= (5,5,0) and symmetry equivalent directions. The rotation direction of the magnetization in 
all magnetic domains, determined by diffraction of circularly polarized light, was found to be unique 
and in agreement with predictions of a combined approach based on a spin-model complemented by 
relativistic density-functional simulations. Imaging of the surface shows that the largest adjacent 
domains display a 120° vortex structure. 

PACS numbers: 75.85.+t, 75.60.Ch, 75.25.-j 



The seminal work of I. Dzyaloshinsky [l| on the rel- 
ativistic origin of weak ferromagnetism in antiferromag- 
netic substances is intimately connected to various emer- 
gent physical phenomena in condensed matter. For 
example, in the skyrmion lattice the very presence of 
antisymmetric exchange interactions (Dzyalonshinskii- 
Moriyaflil) in a non-centrosymmetric crystal stabilizes 
the long period helical structure in zero magnetic field. 
Also, for some spin-driven ferroelectrics (multiferroics), 
the electric polarization is driven by non-collinear mag- 
netic orders; the inverse Dzyalonshinskii-Moriya effect. 
In this case, a phenomenological formulation [3[ shows 
that for cycloidal magnetic structures, i.e. spins rotat- 
ing in a plane that contains the magnetic wave-vector 
(fc), the electric polarization (P) transforms as a prod- 
uct involving the magnetization density and its gradient; 
the so-called Lifshitz invariant of the form P • A, where 
A = (V • L)L — (L • V)-L, and L is the antiferromagnetic 
order-parameter. In a complementary view, the mag- 
netic polarity can be thought of as arising locally from 
spin current Q, as A = k x (Si x Sj), where Si and Sj 
are spins on adjacent sites. Like P, A is a polar vector, 
and will be called magnetic polarity in the remainder. 

In BiFe03, arguably the most studied multiferroic ow- 
ing to room temperature magnetoelectric coupling [EJ, 
ferroelectricity is the consequence of an improper struc- 
tural transition at T c ~ HOOK to the polar space group 
R3c. In bulk samples, the magnetic ordering transition 
occurs at Tat ~ 640K. While the two do not coincide, 
the respective order parameters are coupled through anti- 
symmetric exchange, i.e., P drives the appearance of the 



inhomogeneous magnetization through a coupling term 
7AP, where 7 is a coupling constant, a scenario origi- 
nally proposed by Kadomtseva [fj. The magnetic struc- 
ture can be described locally as canted G-type, but with 
a long period modulation (~ 620 A) in the hexagonal 
basal plane [7j. Subsequent studies j8|, [9[ determined 
that the modulation is cycloidal with the spins rotating 
in the (fc, z)-plane where fc can take the three symmetry- 
equivalent directions fci = ((5,(5,0), k^ = (5,-25,0) and 
&3 = (—25, 5, 0) in the hexagonal setting of the R3c group 
(employed throughout), and 5=0.0045 at 300K. 

In this letter, we study the magnetic domains at the 
surface of a millimeter-size single crystal of BiFeOs with a 
single ferroelectric (FE) domain. Using the high momen- 
tum and spatial resolution of synchrotron X-ray diffrac- 
tion, combined with circular polarization of the beam, 
we determine the absolute rotation direction of the mag- 
netization in individual magnetic domains, which are 
found to have the same magnetic polarity. The sign of 7 
is determined and compared to model-Hamiltonian and 
ah — initio calculations. The large domains observed ap- 
pear to form vortex structures with a closure of the wave- 
vector for three adjacent 120° domains. 

Single crystals of several mm 3 were grown from a 
Bi 2 03/Fe203/B 2 03 flux by slow cooling from 870 °C to 
620 °C. A selected crystal was mechanically cut and pol- 
ished perpendicular to the c-axis, and then annealed to 
remove any induced strain. A piezoresponse force mi- 
croscopy (PFM) of the polished face (not shown) indi- 
cated that the surface had a single FE domain, with the 
electrical polarization pointing down into the sample. We 





FIG. 1. (Color online) Reciprocal-space scans showing mag- 
netic Bragg intensities of the six satellites of the (0,0,9) parent 
reflection, with ki = (5,5,0) (red), ki — (£,—2(5,0) (green) 
and &3 = (—25, 5, 0) (blue) and 5 ~ 0.0045. The x- and y-axis 
are taken respectively along the reciprocal a* direction and 
real space 6-direction. 



label this domain FE|, with the opposite polar domain 
labelled FE^. The synchrotron X-ray experiments were 
performed at Diamond Light Source (UK) on Beamline 
116 [10|. A horizontally polarized beam with a flux of 
~ 10 12 photons per second was delivered by a linear un- 
dulator and tuned to an energy of 5.8 keV, off resonance 
of chemical elements present in BiFeOs. Circular polar- 
ization of the beam was achieved by transmission through 
a 100 jam thick diamond phase-plate, reducing the inci- 
dent flux by ~ 40 %. The diamond crystal was aligned 
to scatter near the (111) reflection in transmission. For 
a certain deviation of AO from the Bragg condition, the 
crystal behaves as a quarter wave plate giving circular 
light. The handedness of the light is determined by the 
sign of A#, which was calculated by dynamical scattering 
theory, and confirmed through experimental calibration 
of the beam line by measuring the X-ray dichroism of a 
standard ferromagnet. 

To prevent contamination of the magnetic signal by 
charge scattering from neighboring structural reflections 
(5 is extremely small), we focussed on magnetic satel- 
lites of the N= (0,0,9) reflection, which is extinct by 
the presence of c-glide planes. Additionally, contami- 
nation from multiple scattering was fully eliminated by 
positioning the sample at an azimuthal angle ^=-170.0° 
with respect to [1,0,0]. Diffraction of A/2 X-rays was 
made negligible by employing up-stream harmonic re- 
jection mirrors. The magnetic signal was clearly iden- 
tified using the full X-ray beam size (100 jam vertical x 
350/im horizontal) with linearly-polarized light scanning 
in reciprocal-space around the positions of the six satel- 
lites N ±k± : N ±k2 : N ±&3, for various positions on the 



FIG. 2. (Color online) a) Photograph of the polished crystal 
surface of BiFeOs normal to the (001) axis (hexagonal set- 
ting, see text for details). A v-shaped defect is seen in the 
center of the surface. The downward direction of the electric 
polarization P determined by PFM is shown (cross) together 
with the reciprocal a*, 6* axis (yellow lines), b) Distribution 
of antiferromagnetic domains with wave- vectors ki = (5, 5, 0) 
(red), k 2 = (5,-25,0) (green) and k 3 = (25,-5,0) (blue). 
The direction of propagation of the cycloidal modulation in 
real-space coordinates for each domain is shown. Each pixel is 
colored according to the diffraction signal (domain) present. 
In the case of multiple diffraction peaks (overlap of domains) , 
the pixels are shaded with mixed colors, respectively. 

crystal surface. The high momentum resolution allows 
the full separation of the six satellites, shown in Fig. [TJ 
in contrast to previous neutron experiments d, [s| . The 
beam size was subsequently reduced using slits to cre- 
ate a footprint of 50x50/im 2 on the crystal surface. An 
image of the magnetic domains (Fig. [2]) was then con- 
structed by step-scanning the sample position with a step 
size of 50± 1/im, recording the intensities of magnetic 
Bragg peaks N + &i,iV + k2,N + k% using rocking-curve 
scans. This procedure lead to the identification of three 
large magnetic /c-domains corresponding to fci, k^ and 
&3, shown in Fig. [2j and to some smaller domains at 
the edges of the scanned surface and around a sizeable 
crystal imperfection in the center of the specimen. The 
three main domains are extremely large, reaching up to 
500 jam in some directions. Note that the average pene- 
tration depth of the X-ray beam is 3.3 jam at this energy, 
placing a lower bound on the domain thickness. Despite 
the long period of the modulation (620 A), this result 
indicates that each domain corresponds to several hun- 
dred magnetic periods. The real space directions of the 
wave- vectors are shown in Fig. [2)3. It appears that the 
modulation of the magnetization follows a 120° vortex 
structure described by the path k\ — > ks — » &2 when 
rotating anticlockwise on the crystal surface. 

To determine the absolute rotation direction of the 
magnetization in each domain (magnetic polarity), scat- 
tering data were collected using circularly polarized light. 
For alternate chiralities of the X-ray beam (left/right 
handed), the intensities of the magnetic signals N + fci, 
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N + &2 5 and TV + fc 3 , were recorded after analysis with 
a pyrolitic graphite crystal as a function of the analyzer 
angle 77, where 77=0 and ?7=90 correspond to the a' and n f 
polarization channels (perpendicular and parallel to the 
scattering plane), respectively. The incident-light polar- 
ization is described by the Stokes vector P 8 =(Pi^P2^Ps) 
PJj, where Pi, P2, P3 represent respectively the degree 
of linear polarization along a and 7r, oblique polarization 
(±45°) and left or right circular polarization. Pi and 
P2 have been determined by fitting the variation with 
77 of the Thomson scattering intensity for the reflection 
(0,0,6), taking into account the cross-channel leakage of 
the analyzer. | ^3 1 was determined by supposing a fully 
polarized beam, i.e. IP3I = y 1 — P x 2 — P 2 2 • In our mea- 
surements, right and left handed light was 93% and 92% 
circularly polarized, respectively (see supplementary in- 
formation for the detailed calculations and conventions 
used). For each magnetic domain, the intensity (Im) of 
the corresponding diffraction peak was evaluated using 
the density-matrix formalism [l2j: 

I M (Q,P S ,V) =tr[D(i 1 ).V m (Q).p(Ps)-V m {Q) ] } (1) 

where p is the density-matrix representing the polariza- 
tion of the incident beam, and D the matrix representing 
the analyser configuration. V m =B.M(Q) is the scat- 
tering amplitude where B is expressed as a two by two 
matrix on the basis of the a and tt polarizations [11| and 
M(Q) the magnetic unit-cell structure factor. For the 
peaks at Q = (0, 0, 9) + ki (z=l,2,3): 



M(Q)=6f(Q)[M\\-piM x ].e 
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where f(Q) is the magnetic form factor for Fe 3+ , calcu- 
lated in the dipolar approximation from [13], M\\ and 
M z are the magnetization vectors of the cycloid along 
ki and the c-axis, respectively, and z is the fractional 
coordinate of Fe in the unit-cell (z=0.2208 at 300K). In 
our conventions /3 = +1 and /3=-l correspond to cycloids 
rotating counterclockwise (CCW) and clockwise (CW), 
respectively, when the structure is viewed propagating 
along ki and c is up. 

Comparison of intensities collected on the three main 
domains and calculations assuming circular cycloids (Fig. 
[3]), unambiguously demonstrates that all magnetic con- 
figurations rotate CW following our definition. This 
is inferred from the 77-positions of the I m extrema ob- 
tained with both light polarizations, which would be in- 
terchanged for a structure of opposite magnetic polarity. 
Within our conventions, A is oriented in the +c direc- 
tion, antiparallel to P. Refining the ellipticity of the 
cycloid (M z / M\\) does not lead to significant improve- 
ments. This, and the failure to observe higher order mag- 
netic satellites, supports the picture of a harmonic mod- 
ulation at 300K, discussed in (13. Il5|. No improvements 
of the fit were obtained by considering a slight tilt of the 
cycloidal plane, as recently suggested [l6| . 
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FIG. 3. (Color online) Top: Variation of the scattered X- 
ray intensity with the analyser angle n (circle symbols) for 
three magnetic reflections (£,£,0), (5 -25,0), (-25,5,0). The 
red (blue) color corresponds to the signal observed with a 
right-handed (left-handed) X-ray incident polarization. The 
solid lines show the results of a least-square refinement of the 
BiFeOs magnetic structure assuming j3=-l (CW, see text for 
details). Bottom: Calculated variation of the scattered X-ray 
intensity with the analyser angle 77 assuming f5—-\-l (CCW, 
see text for details). The direction of electric polarization P 
is shown as a green arrow. 



The relationship between ferroelectric and magnetic 
polarity was further investigated through ab-initio spin- 
constrained calculations in the framework of density- 
functional theory (DFT). The VASP code [17] with the 
PAW pseudopotentials 1 181 was employed within the 
GGA+£7 approach 0, EH (U ranging between 3 and 
7 eV and J=l eV for Fe d-states) including spin-orbit 
coupling, with a plane-wave cutoff of 450 eV. The total 
polarization was calculated via the Berry-phase formal- 
ism [2l|, [22| . Structural parameters for the FE phase were 
taken from Ref. [23j . Due to its long periodicity, the true 
modulation of the magnetization is currently unaccessi- 
ble to DFT. The modulation angle of the antiferromag- 
netic order parameter is given by = 2n(q x x + q y y), 
where q = /ci,&2, or ks. Choosing q = £3, correspond- 
ing to a cycloidal modulation of spins rotating in the 
ac-plane, one needs a supercell na x 2nb x c in order to 
accommodate = 2ir/na. The largest possible super- 
cell, 2a x 4b x c, contains 240 atoms (just within the 
capabilities of state-of-the-art DFT simulations) and has 
modulation angle ir/a. Accordingly, we considered a hy- 
pothetical spin configuration where the cycloidal period 
is reduced to two unit cells along a, with spins rotating 
CW (see Fig. HJ left panels) or CCW. The total energies 
of the two states are then compared in two symmetry- 
equivalent FE states with opposite polarization, and in a 
reference paraelectric (centrosymmetric, R3c) structure. 
As shown in Table U the paraelectric state is degenerate 
with respect to magnetic polarity, which is then lifted in 
both FE states. Furthermore, the energy favored state 
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TABLE I. DFT results obtained for [7=5 eV, J=l eV. The 
energy difference is denned as AE — Ecw — Eccw- FEt and 
FE^ are characterized by opposite collective displacements, 
r, respectively upward and downward, of Bi sublattice with 
respect to layers perpendicular to c axis. 





r (A) 


P c (/zC/cm 2 ) A E (meV/Fe) 


Favored rotation 


FEt 


0.668 


105.17 -2.34 


CW 


PE 










fe; 


-0.668 


-105.17 2.34 


ccw 



switches when polarization is switched. The reliability 
of this trend has been checked for different values of U, 
as well as within a conventional local-density approxima- 
tion, giving \AE\ between 1.1 and 4.7 meV/Fe. These 
findings strongly point to a tight relationship between 
the magnetic polarity of the cycloidal modulation and 
the FE polarization. However, the rather large energy 
difference AE, as well as the disagreement of the pre- 
dicted magnetic polarity with the experimental finding, 
are most probably due to the artificially short modulation 
of the magnetic configuration imposed in DFT calcula- 
tions. Testing this hypothesis by mapping the energy 
evolution as a function of the modulation vector would 
require very demanding - if at all possible - DFT calcula- 
tions. Instead, we adopted a different strategy as follows. 

We introduce a Heisenberg-like spin model with near- 
est neighbor (nn) and next nearest neighbor (nnn) sym- 
metric, as well as antisymmetric exchange interactions. 
The symmetric exchange interactions have been esti- 
mated by mapping the DFT energy of collinear ferro- and 
antiferromagnetic spin configurations onto the Heisen- 
berg model, giving J nnn / 'J nn ~ 0.03, consistent with the 
value extrapolated from spin-wave dispersions 
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The antisymmetric exchange interactions for a given di- 
rection of P are captured through the magnetoelectric 
coupling constants 7 (nn) and 72 (nnn weight). 7 and 72 
are then estimated by imposing the following constraints 
on the mean-field Heisenberg energy: i) the minimum of 
the energy occurs at the experimental modulation angle 
Oexp ~ 3.24° and ii) the energy difference at 6 = tt (i.e. 
the spin configuration simulated in our DFT calculations) 
is equal to AE, as evaluated from first principles. Under 
these assumptions we can estimate 7 2.38 • 10~ 4z V and 
7f £ = 0.6, with \AE(0 exp )\ ~ 0.11 meV/Fe, for U = 5 
eV (the same order of magnitude was obtained for U = 3 
e V and U = 7 e V) . Following Ref . [6| , the inhomogeneous 
magnetoelectric coefficient in the framework of Landau 
theory of phase transitions would be 7 = 4ttA/IP c ~ 
5.8- 10~ 4 V (with exchange stiffness A = 1.87- 10 5 eV/cm 
0, modulation period / = 620 A and assuming the cal- 
culated P c = 105.17 fiC /cm 2 ), in good qualitative agree- 
ment with our estimate. Our model analysis also un- 
derlines the relevant role of nnn interactions, as through 
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FIG. 4. (Color online) Sketch of the considered magnetic 
configuration in the 2a x 4b x c hexagonal cell of BiFeOs. 
Upper left: side view. Bottom left: spin configuration for a 
selected layer of Fe ions. Right panel: mean-field energy as a 
function of the modulation angle for the FE^ domain, with all 
the parameters estimated from DFT with [7=5 eV (see text); 
vertical dotted line marks the experimental 6 e x P , thick (thin) 
lines corresponds to total energy with (without) next-nearest 
neighbor contribution J nnn - A zoom for small modulation 
angles is also shown. The inset shows the energy difference 
between CW and CCW configurations for the optimal 7,72 
and by artificially modifying the nnn contribution 72. 



including J nnn the mean-field Heisenberg energy almost 
reproduces the DFT results even at = n, where the 
only constraint has been imposed on AE (Fig. As an- 
ticipated, the energy-favored magnetic polarity appears 
to depend strongly on the modulation angle of the cy- 
cloidal configuration and on the relative weight of nn and 
nnn antisymmetric exchange interactions, which give rise 
to opposite energy contributions with a different depen- 
dence on (as detailed in the supplementary informa- 
tion). For 7 < 0.7, the energetic competition between 
nn and nnn interactions causes the favored magnetic po- 
larity to change sign when moving from short to long 
modulation periods, therefore reconciling DFT and ex- 
perimental results. 

In summary, magnetic domains of up to 500 /im have 
been observed at the surface of a single crystal of BiFeOs 
consisting of a single ferroelectric domain. The magnetic 
cycloids in each domain were found to propagate with a 
unique rotation direction imposed by the electric polarity 
of the crystal, in agreement with the predictions of our 
theoretical study if nnn interactions are taken into ac- 
count. In future studies, it would be of interest to observe 
the switching of the rotation direction of the magnetic 
cycloids upon switching of the ferroelectric polarization 
by an applied electric field, as observed in TbMnOs [26^ 
and predicted by our calculations. 
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